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Abstract 

Hazard rates play an important role in various areas, e.g., reliabil- 
ity theory, survival analysis, biostatistics, queueing theory and actuarial 
studies. Mixtures of distributions are also of a great preeminence in such 
areas as most populations of components are indeed heterogeneous. In 
this study we present a sufficient condition for mixtures of two elements 
of the same natural exponential family (NEF) to have an increasing haz- 
ard rate. We then apply this condition to some classical NEF's having 
either quadratic, or cubic variance functions (VF) and others as well. A 
particular attention is devoted to the hyperbolic cosine NEF having a 
quadratic VF, the Ressel NEF having a cubic VF and to the Kummer 
distributions of type 2 NEF. The application of such a sufficient condition 
is quite intricate and cumbersome, in particular when applied to the latter 
three NEF's. Various lemmas and propositions are needed then to verify 
this condition for these NEF's. 

Key words: Natural exponential families; mixtures; variance functions; 
quadratic variance functions; cubic variance functions; hyperbolic cosine 
NEF; Ressel NEF; Kummer type 2 NEF. 

1 Introduction 

Hazard rates (also called failure rates) play an important role in various areas, 
e.g., reliability theory, queuing models, survival analysis and actuarial studies. 
Mixtures of distributions are also of a great preeminence in such areas as most 
populations of components are indeed heterogeneous. A comprehensive list of 
references on the behavior of hazard rates for mixtures of distributions can be 
found in the monograph by Shaked and Shanthikumar (2007) and the references 
cited therein and also in Block et al. (2003). 
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In this study wc present (Proposition 2) a sufficient condition for mixtures 
of two elements of tlie same natural exponential family (NEF) to have an in- 
creasing hazard rate. We then apply this condition to some classical absolutely 
continuous NEF's having cither quadratic or cubic variance functions (VF's) 
(c.f., Morris, 1982, Letac and Mora, 1990) and others NEF's as well. A par- 
ticular attention is devoted to the hyperbolic cosine NEF having a quadratic 
VF, the Rcsscl NEF having a cubic VF and the Kummcr distributions of type 
2 NEF. The application of such a sufficient condition can be intricate, in partic- 
ular when applied to the latter three NEF's. Various lemmas and propositions 
arc then needed to verify this condition for such NEF's. Accordingly, we ded- 
icate Sections 4, 5 and 6, respectively, for these three NEF's. In Section 3 we 
consider the rather easy application of the sufficient condition to three NEF's 
having either quadratic or cubic VF's, namely, the normal, gamma and the in- 
verse Gaussian NEF's. Our sufficient condition stems from the following seminal 
result by Glaser (1980): 

Proposition 1 Suppose that the probability density s{x), concentrated on the 
interval (a,oo) (with — oo < a < oo^, is positive such that —b{x) — logs(a;) is 
concave. Then the mapping x i— >■ log s{t)dt is concave on (a, oo) and the 
hazard function h{x) = s{x)/ s{t)dt is increasing. 

(Ron Glaser observes for the one line proof that since b'{x) is nondecreasing, 
one has (l//i)'(x) = e''^^^-^^*'>ib'{t)-b'{x))dt < 0). If the probability density 
is s = e"** the fact that b is convex is by Proposition 1 a sufficient condition 
for having h increasing but not a necessary one: see the remark in Section 2 
introducing the Glaser set as well the Jorgensen and the Karlin sets, or consider 
the density S2 — e^^^ in Section 5 for which h in increasing and &2 is not convex. 

Our sufficient condition for a mixture of two members in the same NEF to 
have an increasing hazard rate is as follows: Suppose that the NEF is written 
as 

{e-^"-'^(^)-'(")l(a,oo)(x)cix, Ae A}, 

where A is an interval and — oo < a. Suppose also that b"{x) > for all x > a 
and denote T{x) = 1/ \Jb"{x). We show in Proposition 2 that if there exists 
c > and S M such that the inequality cT{x) < cosh(c.T -I- d) holds for all 
a; > a, we then can find pairs Ai and A2 = Ai -|- 2c in A and a mixing coefficient 
p e (0, 1) such that the mixture density 

has an increasing hazard rate. This simple condition relics on the fact that the 
mixture is employed with two elements of the same NEF. However, this two 
element mixture result is apparently not extendable to a more multi-element 
mixture situation. 
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2 A sufficient condition for mixtures of mem- 
bers of the same NEF to have an increasing 
hazard rate 

Consider an absolutely continuous NEF concentrated on (a, oo) with — oo < a < 
oo, and generated by a locally integrable function s on (a, oo) such that 

L{X) ^ e'^f^) = / e-^^s{x)dx, (1) 



the Laplace transform (LT) of s{x), exists on a nonempty open interval A. The 
corresponding NEF is then given by the set of probability densities on (a,oo) 
of the form 

{exp{-Aa;- fc(A)}s(a;)da;, A e A}. (2) 
Let z/((iA) be a probability on A and suppose that the function on (a, oo) defined 

exists. Thus, s{x)R{x)dx is a probability density on {a,oo) and it is a mixture 
of the elements of the NEF. This probability density has the hazard rate 

_ six)R{x) 

L s{t)R{t)dt 

Proposition 1 shows that hu{x) is increasing ii x \og{s{x)R{x)) is concave, 
or, equivalently, if s{x) > for all a; > a, if s"{x) exists and if on (a, oo) one has 

s"{x)six)-{s'ix))^ R"{x)R{x) ~ {R'jx))' 

^) + m^) - 

We now have the following proposition when is a mixture of two Dirac mea- 
sures. 

Proposition 2 Consider the special case of the hazard rate in ^ with 

= pSxi + {I - p)S\2, (6) 

where p £ (0,1) and Ai < A2 with Ai and A2 in A. Assume that on (a, 00), 
s{x) > 0, —b{x) = logs(x) is concave and that s"{x) exists and define 

T{x) = (7) 

= pe-MAi)^ P2 - (1 - p)e-^(^^), c = and d - log VpIM- (8) 

Then the hazard rate with v as in (0) is increasing if for all x > a 

cT(x) < cosh(cx + d). (9) 



3 



Proof. The proof is a straightforward apphcation of ([S]). Indeed, for in ([5]) 
and R defined by ([3]), 

i?(a;) =pe-*-'(^i'e-^i^ + (1 p)e-'''(^^)e-^^^ = pie'^^^^ + ^26"^=^=", 

implying that R"{x)R{x) - ^ =pip2(A2 - A^^e^^^i+^^js:^ Accordingly, 

the inequality ([5]) becomes for this particular case 



P1P2 



(A2 - Ai)2e-(^i+^^)- < (pie-^i- + P2e-^n^ (-'^^] ' , x > a. (10) 



Since logs(a:) is concave, {—s'{x)/s{x)y > on (a, 00), so that T{x) in ([7]) is 
well defined. Thus, with the notations (|S|), the inequality ((TU]) is equivalent to 

2cT{x) = T(x)(A2 - Ai) < . fEle'^^^ + 1?!^-^^^ = 2 cosh(ca; + d). 

V P2 V Pi 

which is dni). ■ 

In the next sections, we are going to consider a number of absolutely con- 
tinuous NEF's on the real line, generated by a density s, and check for each of 
them whether T exists or not, that is whether s is log concave or not. When T 
exists we will have to discover which (c, d) with c > are such that ^ holds for 
all X. As we shall see, for some NEF's such that s is log concave it may occur 
that ([9]) does not hold for any (c, d). 

The system ([8|) of equalities links the three parameters (Ai,A2,p) with the 
two parameters (c, d). Suppose that we are given a pair (c, d) satisfying we 
therefore may choose arbitrarily the mean A = ^^^^^ in A such that Ai = A — c 
and A2 = A + c are in A. Having done this choice of A, the value of of the mixing 
coefficient p in ([8|) can be determined exactly as 

^ e'^L{\-c) + e~'^L{X + cy ^ ' 

where L is the Laplace transform (LT) of the generating density s (recall that 
s is not necessarily a probability). We note, however, that the LT is not always 
expressible in terms of simple functions but rather in terms of transcendental or 
implicit functions, in which case a numerical search is then needed to find the 
(c, d) interval on which the appropriate mixture density possesses an increasing 
hazard rate. As this paper is rather theoretical, we do not intend to pursue such 
a numerical search. 



Remarks on the Jorgensen, Karlin and Glaser sets. Given a density s 
on (a, 00) with LT ^ such that A is not empty, the set J{s) of a > such 
that is still a LT of a positive measure jia is called the Jorgensen set of s 
(see for instance Letac and Mora (1990) and the references cited therein). By 
definition, J(s) is a closed additive semigroup. Note that s generates an NEF 
of infinitely divisible distributions if and only if J(s) = [0, 00). If not, J(s) can 
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be complicated. For instance, a consequence of the short and elegant paper by 
Ben Salah and Masmoudi (2010) is that J(s) — [l,oo) if 

1 3 

s{x) = -e""^l(o,oc)(a;) + -e~''+^l(i^oo)(a;). 

Note also that fia could have a continuous singular part for some small a S 
J(s) although appropriate examples are rather complicated. If J* (s) C J(s) 
is the set of a such that fj,a(dx) = Sa{x)dx has density we trivially have that 
J*(s) + J(s) C J*{s), since the convolution of a measure with a density with 
any measure has a density. 

Consider now a non trivial result due to Karlin and Proschan (1960) (see also 
Karlin, 1968, p. 152 and Barlow and Proshan, 1965, p. 100) which says that if s 
and i are probability densities with increasing hazard rate, then the convolution 
s * ^ has the same property. Therefore let us introduce the Karlin set K{s) of 
a G J*{s) such that Sa has an increasing hazard rate. The above property 
shows that K{s) is a closed additive subsemigroup of J*{s). For instance if 
s{x) — e^^l(o,oo)(a;), it is a simple exercise to see that K{s) = [l,oo). 

Finally, consider the Glaser set G(s) of s which is the set of a in the Karlin 
set K{s) for which the conditions of Proposition 1 are met, that is such that if 
Sa = e~^°' then ba is convex. Although in many cases G{s) coincides with K{s) 
the Glaser set G{s) is not a semigroup. Indeed, we shall face with an example in 
Section 4 relating to the Ressel distribution si where G{s) is a bounded interval; 
thus distinct from the semi group K{s). Usually the Karlin set is more difficult 
to find than the Glaser set. 

3 Applications related to NEF's with quadratic 
or cubic VF's (normal, gamma and inverse 
Gaussian NEF's) 

As already noted in the introduction, quadratic VF's include six NEF's of which 
only three have densities: Normal, gamma and hyperbolic cosine (c.f. Morris, 
1982). Cubic VF's include also six NEF's of which only two have densities: In- 
verse Gaussian and Ressel (c.f. Letac and Mora, 1990). The set of our examples 
will include all of the five absolutely continuous NEF's having either quadratic 
or cubic VF's and also another one, the Kummer distributions of type 2 NEF. 
The present section deals with the normal, gamma and inverse Gaussian NEF's 
and Sections 4, 5 and 6 consider the three other ones. In what follows and 
whenever feasible, we provide, for each of the examples, with their respective 
VF (F, ri), where V is the VF corresponding to ([2]) and is the domain of 
means. 

Example 1: The normal NEF 

The normal NEF has a constant VF, i.e. {V,Vi) = (cr^M). For a fixed 
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standard deviation cr, the generating density is 

1 _ x2 

s(x) — , e 2^, a = — oo. 



Trivially here the Glaser set G{s) is (0, oo). This leads to T{x) = a and fc(A) = 
2^^. The inequality (|9]) is fulfilled for any a; G M if and only if cu < 1, or 
equivalently, if |Ai — A2I < 2/ct, a result that was already obtained by Block et 
al. (2005). 

Example 2: The gamma NEF 

The gamma NEF, concentrated on (0,oo), has a VF y(/i) — a~^fi'^ and 
fl = R"*", where a and /it are, respectively, the shape and mean parameters. For 
a fixed shape parameter a > 0, the generating measure is 

Sa{x)dx = ——1^0 oo){x)dx. 
i (a) 

We now consider three exhausted cases relating to the values of the parameter 
a: a = 1 (the exponential case), a < 1 and a > 1. These observations imply 
that the Glaser set G{si) is [l,oo). 

1. a = 1. Here, s{x) = 1 so that for any v the function log R is convex and the 
inequality ^ (as well as (O) cannot be fulfilled unless v is concentrated 
on one point. 

2. a < 1. Since for this case both log\R and logs are convex, the inequalities 
© or dH) cannot be fulfilled. 



3. a > 1. Here, T{x) = x/\/a ~ 1 and we have the following proposition. 
Proposition 3 For a > 1, the probability density 

/(^) = ^x"-i(pA?e-^^- + (1 - p)X^e-^^n, 2: > 0, 
T(a) 

where Ai < A2, has an increasing hazard rate h{x) — f{x)/ f{t)dt if 

l/a 

Al ~ Vl -P/ 

where 



do = log - — — 

Va - 1 (1 + V") 

Example. For a = 2 the result specializes to the following: since do — log 2 — 2 
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and since 6^^0/2 = e/\/2 = 1.92... we can claim that the following mixing of two 
gamma densities 

f{x) = pXjxe-^''' + (1 - p)Xlxer^^'= 

where Ai < A2 has an increasing hazard rate if (^^^ < i3];(l-92..). For 

instance choosing A2 ~ 2Ai imposes a heavy weight p on Ai, namely .675 < p < 
1. 

Proof. Since c > for studying the inequality ([9|) we write t = cx + d. Thus 



^ becomes: for all t > —d, 

t-d 



< cosht. (12) 



For a fixed a > 1, we determine the set of d values such that (fT2|) holds. Since 
1 1— >■ cosht is a convex function, we look for the point (to, cosh to) such that the 
tangent to the curve cosh has slope l/^/a — 1. Thus, 

1 , l + , , 2a 

smh to = ^^^= , to — log ^^^= and cosh to 



^/^^' ^V^^ (l + V^)^/^^■ 



The equation of this tangent is y ^ {t — do) /y/a — 1, where do is such that this 
line goes through the point (to, cosh to). This implies that 

1 + ^/a 2a 



do = to- (coshto)\/a - 1 = log ~ — — (13) 

y/a-1 (1 + V") 



Such results show that t/\/a — 1 < cosh(t + d) for all t > if and only if 
d > do- The application of this fact is that (fTO|) holds for all a; > if and only 



Ai / 1 - p^ 



if \fp\fp2 > e"*", or, equivalently, if 



A 
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Example 3: The inverse Gaussian NEF The inverse Gaussian NEF has a 
VF V{ii) = a^^^i^ with VL — M"'', where a > and n is the mean parameter. 
Here a = and for a fixed a > 0, the corresponding NEF is generated by 

six) = — ^x^^/^e"'^. 
V2^ 



This implies that {— s' {x) / s{x))' = (3a; — 2a^)/2a;^, which is not a positive 
function. Thus Proposition 2 is not applicable. 
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4 The hyperbolic cosine NEF 



The hyperbolic cosine NEF Ta has a VF y(^) = ^ + a with 17 = R, where 
a > 0. The convex support of Ta is K (i.e., a — — oo). The generating measure 
of J- a is 



^ 2 ' 



^ dx 



r(a)- 

(See Morris 1982, for details). Moreover, its LT is defined on A = (-f , |) by 
.^q(A) = (cosA)~" (since it is 1 for A = 0, this shows the non obvious fact that 
fia is a probability) . The hyperbolic cosine fii distribution is more known as the 
hyperbolic secant (hereafter, HS distribution or HS NEF). Various probabilistic 
properties of the HS distribution have been derived, though it is rarely used 
in applied statistics, probably due to its intricate structure. Although this 
distribution is not used much in applications, it does, however, have two curious 
features: Like the normal distribution, the density of /ii is proportional to 
its characteristic function; the sample mean and median are, asymptotically, 
equally efficient. A probabilistic interpretation of fix is available: consider a 
standard complex Brownian motion Z = X + iY with Z{0) =0 and the hitting 
time T of the set {x + iy ; \y\ > 7r/2}. Then X{T) ~ /ii : to see this, consider 
the process M{t) = expsZ(t). Since z i— )■ e^^ is analytic, it is harmonic, M 
is a martingale and E(M(r)) — 1 gives the desired result. A newsworthy 
statistical analysis and data fitness can be found in Smyth (1994) and recently 
in Sibuya (2006) (a complete English version of the latter paper is available by 
corresponding the author). 

Denote by Sa — e^^" the density of /Iq,. The fact that the function ba is 
convex if and only if a > 1 has been proved by Shanbhag (1979). We give a 
different proof in the following proposition: 

Proposition 4 The function ba is convex if and only if a > 1 (in other terms 
the Glaser set G'(si) is [l,oo)). More specifically for a > \ we have 

{a -I) e^"(cosu)"^2 ^u. (14) 



1 



Sa{x) 



-Tr/2 



and for a < 1 the function b'^ is negative in the interval 




Proof. Formula ((141) is the particular case v = — 1 > and a = ix the 
classical formula 

{cosuY cos audu= , xi r (15) 

which can be found in Gradshteyn and Ryzhik (1980), page 372, 3.631 formula 
9. Now (jl4p shows that 1/sq is the Laplace transform of the positive measure 

[a - 1)(cosu)"~21(_^/2,^/2)(m) du (16) 
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which implies that the function ba = —logs a is strictly convex for a > 1. 
For a = 1 we can see that bi — — logsi is convex by the same trick since 
1/si — 2 cosh ^ is the Laplace transform of the positive measure 5_7r/2 + 
Note that is the weak limit of (IT51) when a 1. 

Suppose now that < a < 1. We use the digamma function V = T'/F and 
its derivative. If z is a complex number with positive real part : 



^ {n + z) 



A easy calculation leads to 



where for t > we define 



C — t v-^ 
mJt) = ^ + > 

(c' + tf ,tt((n + c)2 + ir 

For showing that for c < 1/2 the function 1 1— >■ (pc{t) is negative on some interval 
observe that for fixed t > the function on (0, oo) defined by w i— ?> (u+t)'^ 
decreasing when u > 3t. As a consequence if (1 + c)^ > 3t wc can write the 
majorization of the sum of a series by an integral 



[c^ + ty Jo {{v + cy + ty (c^ + ty 

(here we have used (-j2^/-)2 = ~'^~^^t^)- This shows that ^cif) < when 

£!il±£) < t < ii±£l^. Since c < 1/2 we haveii±£l^ - %t£) = _^(i_4c2) > o 

and this interval is not empty. Replacing c by a/2 and t by a;^/4 we get that 
is negative in the interval indicated in the statement of the proposition. 
■ 

For part < a < 1 our proof of Proposition 4 is elementary. For part a > 1 
our proof is based on the formula [15] and Laplace transforms. The compact 
and ingenious Shanbhag's proof rather relies on Fourier transforms through the 
formula b'/^{x) — e^*^^'^da{t)dt where 

This formula is derived from an integral formula for ba which can be found in 
Zolotarev (1967) and which is obtained from the Levy measure of the infinitely 
divisible distribution of log X when X is "fa/2 distributed. If a > 1 the function 
da is an integrable characteristic function (corresponding to a Cauchy distri- 
bution with parameter {a — l)/2 convoluted with the density 7r/2(cosh7ra;)^)). 
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The Fourier inversion formula shows that > for all a:. If < a < 1 we 

have da(x) > 1 around zero and da cannot be a characteristic function. This 
prevents 5^ to be positive by a careful but standard reasoning using again the 
Fourier inversion and this concludes the Shanbhag's proof. 

Since the functions ha and T = 1/ yj^l are not simple when a is not an 
integer we therefore emphasize the analysis of the respective mixtures for the 
two cases: a = \ and a = 2. In principle an analysis similar to the case a = 2 
below could be also performed for a = 3, 4, . . . but the case a = 2 is creative 
enough to let us think that higher cases are difHcult. 

The hyperbolic case a = 1: The most popular member of the J-"q's is re- 
lated to this case. As mentioned in the proof of Proposition 4 the corresponding 
density for a = 1 is 

^^^") = 2jh^- 
Thus the above results are applicable to this J^i . More specifically 



si{x) J 2 cosh" 



and thus T{x) = ^ cosh^. In order to study the inequality ([9]) for this partic- 
ular case, we use the following lemma. 

Lemma 5 Let a and u be positive numbers and v be a real number. Then, the 
following inequality 

a cosh a; < cosh(ua; -(- v) 
holds for all real x if and only if a ^ (0, 1], w > 1 and \v\ < Vq = Vo{a,u), where 

vo = ulog (- + ) - log {A + B) 



with the notation A = \ °\ ci'^d B 



Proof. . Letting x 00, we have that a cosh a; ^ ae^ and cosh{ux+v) ^ e"^, 

implying that u> I. Letting x = —v/u shows that acosh(— u/w) < 1 and thus 
o < 1. In the sequel we assume that m > 1 and treat the case u = 1 separately 
after. Now, we introduce the two positive numbers xq and vq such that the 
two curves x 1— > a cosh a; and x 1— > cosh('ua; — vq) are tangent on a point of the 
abscissa Vq- Thus, they satisfy the two equations: 

acoshx — cosh(MXo — vq) and asinhx — wsinh(ua;o — vq). 

Squaring these two equations and using the fact that cosh^ t — sinh^ ^ = 1, we 
get a linear system in cosh^ xq and cosh^(MXo — vq), whose solution is 

cosh^ xo — and cosh^fuxo — vq) = A^. 
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Since t > and y = cosht, it follows that t = log(y + yjy'^ — I. Thus, xo = 
log (4 + ^) and Wo = ulog(4 + ^) - log (A + S) (note that mxq - wq > 0). 
To end the proof of we show that a cosh x < cosh(Mx + v) for all real x 
would imply that \v\ < vq. Since the function v i— >■ cosh(ua;o — v) ~ acoshxo is 
decreasing on the interval (—oo,uxo) and is zero on vq (which belongs to this 
interval), we get that v < vq when cosh(ua;o — v) — acoshxo > 0. Similarly, 
because of the symmetry of cosht we have —vq < v. 

We now prove . Assume that a < 1 < u and \v\ < vo- Denoting f{x) = 
cosh(wa; — v) — a cosh a;, then since f{xo) = /'(xo) = 0, the Taylor formula gives 

fix)= r{x-t)f"(t)dt. (18) 

J Xo 

We use the latter formula to show that f{x) > for x > xq. Note that since 
f"{x) > f{x), then f"{xo) > and ([T^ implies that f{x) > on some interval 
(a;o,a:i). Now suppose that there exists X2 > xq such that f{x2) = 0. Without 
loss of generality we assume that f{x) > on {xo,X2)- Thus f"{x) > on 
{xo,X2)- Since by (fT8|) . f{x2) = is impossible, we obtain that f{x) > for all 
X > xq. To prove that f{x) > for all x < a;o is similar. 

We now consider the particular case u = 1. The inequality a cosh a; < 
cosh(a; + w) is equivalent to e^^(e" — a) > a — e^". By letting x — ±oo, it can be 
easily seen that the latter inequality holds for all x if and only if |w| < — log a. 



We do not apply the full strength of this lemma for our problem, but instead 
study when the inequality 

2c , Tra; , , ,. 

— cosh — < cosh(ca; + a), 

TT 2 

holds for all x. 

To fit with the notations of the latter lemma, denote t = ^ which leads to 
— cosht < cosh(— a; + d). Now, the lemma implies that if this inequality holds 
for all t then a =^ < l< u = ^. Thus — = 1, but since c = i(A2 - Ai) we 
must have A2 = Ai + tt. However, this is impossible since the corresponding LT 
Li(A) = (cosA)^^ is not defined outside of the interval (— ^, |). To end up this 
discussion, no mixing can give increasing hazard rate for the NEF generated by 
the density si{x) = ^ ^ . 

The hyperbolic case a — 2: As we are going to see here, the situation is 
more favorable when dealing with the following direct consequence of (|14p : 

S2{x) 



2 sinh ^ 

Proposition 4 has shown that 62 = — log S2 is a convex function, so we are in 
position to use Proposition 1. We have the explicit calculation 



sinh^ ^ 
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Since |t| < |sinht|, clearly b2{x) > and we have a direct proof of the log 
concavity of S2- Thus we have to study the set of (c, (i)'s such that the inequality 

X sinh 

cT{x) = r 2 < cosh(ca; + d) (19) 

y'sinh^ ^ - {^y 

holds for all real x. For this we use the following lemma. 



Lemma 6 For all real t we have 

t sinh t 



\J sinh^ t — t'^ 
where an equality occurs when t = 0. 



< V3 + <2 



Proof. The proof follows from the inequality sinh f — — y > 0, which is 
deduced from the expansion of 

sinh^ t-f- — = - cosh2t ---t'^- — > 0. 

3 2 2 3 ^ (271 ! - 

n— 3 



In terms of the function T, the lemma is equivalent to stating that for all x we 
have 

2 



T(.)<-^3+(-)^ 

Since it is difficult to find all (c, d) such that ([9]) holds we shall content to study 
the set of (c, d) such that c|-i/3 + (^)^ < cosh(ca: + d) holds for all a:, or 
equivalently, by introducing k = 2c/ tt and u — + d, to study the set of 
{k,d) such that 

A/3fc^ + (u — d)^ < cosh M 

holds for all u. 



Lemma 7 For fc > and d real y^3k^ + {u — d)^ < coshu holds for all u if 
and only if \d\ < do = \/2 — log(l + \/2) = 0.532... and 

3fc^ < (2 — cosh^ Ud) cosh^ Ud, 

where Ud is the solution of the equation sinh2w = 2(u — d). In particular, the 
inequality \/3fc^ + < cosh u holds for all u if and only if k < y | . 



Proof. The inequality -\/3fc2 + (u — dp < coshu implies |it — dj < coshu for 
all u. Now the minimum do of coshw — u is attained at log(l + a/2) which is 
the solution of the equation sinhu —1 = 0, and thus, do = \f2 — log(l + \/2). 
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Similarly, the minimum of coshu + it is attained at — log(l + and is do- 
Since, — coshu — u < —do < dci < coshu — u, we get that |u — dj < coshzt 
for all u if and only if |d| < do. Now, fix d e [—do, do], then for finding all 
k such that 3fc^ + {u — d)^ < cosh^ u, we look for the (positive) minimum of 
u I—)- cosh^ u — {u — dY ^ which is attained at the point Ud- Letting d = gives 
= and also entails the final result. ■ 

Practical conclusion for a — 2. The Laplace transform of S2 is 1/ cos^ A for 
A e (— f , f )■ According to Lemma 7, we fix any d such that |d| < 0.532.. and 

c> such that c < f = 1.282.. We now choose an arbitrary number A such 
that 

--<Ai = A- c<A + c=A2<- 

and we use formula (|lip for defining the mixing coefficient p depending on d, 
Ai and A2. With this choice, the density 



2sinh^ 



[pe-^i^^ cos^ Ai + (1 - p)e-'^^ cos" A2] 



has an increasing hasard rate. 



5 The Ressel NEF 

Consider the density Sq on the positive real line defined by 

s^{x) - e-''"^-) = ^ ' . , X > 0, a > 0, (20) 
r(.T + a + 1) 

then Sq is called the Ressel or the Kendall-Ressel density with parameter a. It 
is infinitely divisible (in other terms the Jorgensen set J(si) is [0, cxd)) and 

This density appears in various areas. For an M/G/1 queueing system with 
arrival rate A, it is the limiting distribution, as A — > 00, of the length of the busy 
period T{a) ~ a initiated by the virtual time quantity a > (c.f. Prabhu,1965, 
pages 73 and 237). In their characterization of the regression of the sample 
variance on the sample mean, Fosam and Shanbhag (1997) showed that such a 
regression is cubic on the sample mean for only for six distributions, of which 
one is the Kendall-Ressel distribution. Kokonendji (2001) also revealed this 
distribution on his investigation of first passage times on and 1 of some Levy 
processes for NEF's. Additional references can be provided here regarding the 
Kendall-Ressel distribution like Fakes (1995) formula (4.1), but one of the most 
detailed reference is Letac and Mora (1990) who characterized all NEF's having 
cubic VF's, of which, of course, the Ressel NEF is one of them. 

The Ressel NEF generated by ^ has a VF (F, VL) = (^t^{l + ^), (0, 00)^ . 
We are interested in the values of a such that h'^{x) > for all x > 0. One 
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can consult Proposition 5.5 of Letac and Mora (1990) for checking the puzzling 
formula Sa{x)dx = 1 and page 36 of this reference for learning why this 
density can also be called the Kendall-Ressel density. 

Proposition 8 Let 

a—1 (2 — a)x — a"^ + a x (a'^ — I) + (a — a^)x + (2 — a)x'^ 

x'^ x{x + a) X (x + a + l) 

(21) 

then for all x > Q we have h{x) < b'^{x) < g{x). Furthermore there exists 
a number a £ (1.77,1.91) such that b'^{x) > for all x > if and only if 
a € [l,a]. In other terms the Glaser set G(si) is [l,a]. 

Proof. For x > 0, we use the digamma function -0 — F'/F and formula (HZ]). 
Using this notation we have 

f{x) = bl{x) = -^—^ + ^^ + ^'{x + a + l) 



a—1 a 



1 °° 1 

- + ^7 (22) 



X x^ ^ (n + X + a — 1)^ 

n— 2 ^ ' 



Now observe that 



1 _ r°° dt ^ 1 1 r°° dt 

1 + x + a J2 {t + x + a-1)'^ ^(n + j: + a-l)2 x + a~ {t + x + a 

This gives the desired inequalities h < f < g. Clearly the function / is positive 
if a = 1. If a < 1, the function / is equivalent in a neighborhood of a; = to 
which tends to —00. It is obvious that for a > 2, g{x) becomes negative 
ultimately (if a = 2 then g{x) = {2 — x)/x'^{2 + x)). Hence assume that a < 2 
and let us study the sign of g. Since g{x) =0 if 

(a^ - a) + (-1 + 2a- a^)x + (2 - a)x^ = 0, 

then this equation has at least one solution if 

D{a) = (-1 + 2a + a^f - 4(2 - a){a^ - a) = I + Aa - Qa^ + a"" 

is nonnegative, which is the case for a > a* — 1.90321. One of the two possible 
solutions of g{x) ~Q is 



xi = : — > — — > 0, a < a < 1. 

A -2a A -2a 

So .9(2:1) — and hence g{x) < for some a; > 0. On the other hand, we use 
inequality h < f and study the sign of h. Let 1 < a < 2, then if D{a) = 
{a — 1)(— 8 — Aa + 3a^ + a^) < 0, h{x) has no zeros at all and hence h{x) > 
for all a; > 0. Now, D{a) = if a G {-3.48929,-1.28917,1.0,1.77846} and 
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D{a) — >^ oo as a — > oo, hence D{a) < (i.e., / is positive for all a; > 0) if 
1< a < 1.77846. ■ 



Remarks. For studying the log concavity of the density Sq, of the Ressel distri- 
bution one can be tempted to imitate Proposition 4 and to wonder if for a > 1 
the function 

— -- = :reT(x + a + 1) X x'"" 

defined on (0, oo) is the LT of a positive measure. Proposition 8 has shown 
that is is impossible if a > a. This can be explained by the fact that the factor 
X H-> x~^ is the only factor in 1/sa which is not a LT (this observation follows 
from the fact that x is the LT of a stable law with parameter 1, and the 

reciprocal of the LT of a non Dirac measure cannot be a LT). In terms of the 
Glaser and Karlin sets, the density si is quite interesting. Proposition 8 has 
shown that G{si) = [1, a] C K{si) and a striking consequence is that while the 
density si is log concave the density S2 — si * si is not. This demonstrates 
the difference between the Glaser and the Karlin sets. The additive semigroup 
generated by [1, a] is [1, a] U [2, oo) C K{si). One can reasonably conjecture that 
K{si) = [l,oo). 

From Proposition 8 it follows that if a e [1, a], we are led to consider T{x) = 
1/ -^/b" (x) and study the set of the pairs (c, d) such that cT{x) < cosh(ca; + d). 
Since h{x) < b'^{x) = —{s'^/sa)'{x), a sufficient condition for these {c,d) values 
is that cj ^h{x) < cosh(cx + d). Even this simplified inequality is still too 
complicated and we shall content here to consider only the case a = I. For 
this case, we search for the set of (c, (i)'s such that for all x > we have 
c\/x + 2 < cosh(cx + d). Therefore, the next proposition is devoted to the case 
of the NEF generated by the probability density 

si{x) = Y^^^^^{o,oo){x)dx. (23) 

Proposition 9 For c > and d real consider the function 

ip{x) = 4j cosh^(cx + d) — X — 2 

and define xq = log(c + Vc^ + 1) — d) . Then ip{x) > for all x > if and 

only if 



1. either xq < and coshd > y/2c; 

2. or Xo > 0, c < and xq < ^{1 + Vc^ + 1) - j- 



Proof. Since (p'{x) ~ i sinh2(cx + (i) — 1, (p'{xo) < for x < xq and (p'{xf)) > 
for X > Xo, it follows that for xo < 0, (p{x) > for x > if and only if (f{0) > 
or if and only if coshd > y/2c (recall that xo < implies d > 0). Similarly 
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for xo > 0, we have (p{x) > for a: > if and only if if{xo) > 0. Since 
^ cosh^(cj:o + d) — 2c^(l + \/c^ + 1), the inequality ip{xo) > is equivalent to 
Xo < ~^ Vc^ + 1) — J, which can be realized only if the right hand side is 

nonnegative, that is if c < ■ 

Now, having pairs (c, d) to our disposal, we have to compute p as given by 
the formula ([TT]) and we need for this the values of the Laplace transform at 
points Ai = A — c and A2 = A + c. For the Ressel distribution si, its LT 

/•oo 

Li(A) = / e-^''si{x)dx (24) 
Jo 

cannot be expressed explicitly. However, a numerical or graphical calculation 
of Li{X) for a given positive value of A is easily done by means of the following 
proposition. Its statement is equivalent to formula (11) in Fosam and Shanbhag 
(1997) which relies on Prabhu (1965, p. 73 and p. 237). We give an independent 
proof here for sake of completeness. 

Proposition 10 Consider the bijection f from [1, 00) to itself defined by f{x) = 
X — log a;. Then for a given A > 0, the number l/Li{X) defined by \24-^ satisfies 
/(l/Li(A)) = l + A. 

Proof. Let {Y{t))t>o be the Levy process such that for t, A > we have 
£(6""^^^*)) = (1 + A)~* (such a process is usually called the gamma process). 
Define T = min{t : t — Y{t) = 1}. The random variable T— 1 has the Ressel dis- 
tribution (|23l) ( see also Letac and Mora, 1990, p. 27 for this observation). Fur- 
thermore, Theorem 5.3 in Letac and Mora (1990) states that if £(6*^^) = e-«(^) 
andE(e^(i-'^(i)) = e-'^C') then r{q{e)) = 9. Thus, since Li(A) = E(e-^(^-i)) = 
gA-g(-A) and since e-''(®) = we can write r((7(6')) ^ \og{l + q{e)) - q{e) = 9. 
This leads to 

log(l + qi-X)) = qi-X) - A = - logLi(A) 

The elimination of q{—X) between these two equalities gives jj^^ +logii(A) = 
1-f A. Since, for A > we have Li{X) < 1, we get /(l/Li(A)) = 1 + A. ■ 

Practical conclusion for a — 1. We fix a pair (c, d) such that either condition 
1 or condition 2 of Proposition 9 holds. We choose a number A such that 
0<Ai=A — c<A-|-c=A2. We compute numerically Li(Ai) and Li(A2) with 
Proposition 10. The mixing coefficient p is therefore determined by pT|) . The 
density on (0, cx)) 



r(a; + 2) 

has an increasing hazard rate. 



+ (1-P)- 



Li(Ai) ^ "'Li(A: 
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6 The Kummer type 2 NEF 



Let a, A > and b real and consider the number 

fco a-l 

As a function of A, C is proportional to what is sometimes called in the litera- 
ture the confluent hypergeometric function of the second kind or a Whittaker 
function. If 



then the probability s{x)dx = K^^\a,b,\){dx) is called the Kummer distribu- 
tion of type 2 with parameters (a, &, A). Needless to say, that if (a, h) are fixed 
the model {s{x)dx, A > 0} is an NEF. If 6 > 0, this model is generated by the 
beta distribution of type 2, i.e., by 

1 x'^~^ 

Kummer distributions have been studied by Ng and Kotz (1995). Statistical 
aspects of Kummer distributions for waiting times and exceedance statistics 
have been considered by Fitzgerald (2002). The Kummer distributions of type 
one belong to NEF's generated by the ordinary beta distributions. Since they 
are concentrated on the bounded set (0, 1), they are not relevant for our study. 
Accordingly, we study the NEF generated the Kummer distribution of type 2. 
Its VF cannot be expressed explicitly. However, the important fact about such 
an NEF is the formula (P7|) below which gives the LT of /3*^^^ (a, h) in terms of 
the confluent hypergeometric function defined for real a and h such that b is not 
in the set — N. This LT is then given in terms of the entire function 

n— 

Here, (a)o = 1 and (a)„+i = (a + n)(a)„. This formula states that if a > and 
h is not in the set Z of relative integers, we have 

C(a, b, A) = ^^^^ iFi(a ; 1 - 6 ; A) + r(-6)A'' i^i(a + 6 ; 1 + 6 ; A). (27) 
r(a + b) 

In (j27p . the mapping z i-^ l/r(z) is an entire function which coincides with the 
ordinary 1/T{z), 2; > 0, and (|27p can be extended to the case where & G Z by a 
limiting process. The identity (j27p is by no means elementary and its proof by 
the Barnes formula can be found for instance in Slater (1960), formula 3.1.19. 
A probabihstic proof would be desirable. 

To exemplify the use of (|27)) . observe that if a, fe, A > and \iX ^ 7(6, A), Y ~ 
K^'^\a, b, A) and Z 7(0 + b, A) are independent, then X + Y and Z/{1 + Y) 
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have the same distribution K^'^\a + b, —b, A). In order to prove this, for suitable 
t's just consider the LT E(e-*('^+'^)) and the Mellin transform E{Z*/{1 + F)*), 
and then use ([?7)) . 

Proposition 11 If s is defined by S25\) then — (s'/s)' > for all x > if and 
only if I < a and 6 < — 1 with a — 6 — 2 7^ 0. 

Proof. If ^ = a — 1 and B = — 6 — 1, we get 

Trivially, -(s'/s)' > if A > and B > with A + B ^ 0. Conversely, if 
A + 2 Ax + Bx'^ > for all a; > 0, then letting x 00 shows that B > 0. 
Also, letting x —?' shows that A > 0, while A = B = Q would imply that 
-(s'/s)' = O.Therefore if ^ = a - 1 > and B = -5 - 1 with AB ^ 0, we are 

allowed, by using p8|. to consider for x > 0, 

T{x) = l/\/-(-)'(x) " + 



s y/A + 2Ax + Bx2 

We then have to investigate which numbers c > and d real are such that for 
all X > 0, we have 



X + x"^ 



\J A + 2Ax + Bx' 



< cosh(ca; + d) 



For simplicity, we are going to treat only the particular case A — 0, and therefore 
to study the NEF 

I (I + -Ax 

where i? is a fixed positive constant. For this particular case we are looking for 
the values of (c, d) with c > such that for all x > we have cosh(cx + d) — 
X - 1 > 0. ■ 

Proposition 12 For B,c > and d real, consider the function defined on R 
by 

(p{x) = cosh(ca: + d) — x — 1 

c 

and define xq ^ (log ^^^j^^^ ^ . Then (p{x) > for all x > if and only 

1. either xq < 0, a/B < c and coshd > 

2. or xo > 0, c < \/B + 1 and xq < I - 
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Proof. We study the function tp in an elementary way: We get that ^p'{x) = 
y/B sinh(ca; + d) — 1 satisfies ^p'{x) < Q ioi x < xq and ^'{x) > for x > xq. If 
xo < then ip{x) > for ah x > if and only (p{0) = coshd — 1 > and 
this proves part 1. If xq > 0, then if{x) > for all x > if and only 



(^(xq) = I- Xq>0 

c 

which proves part 2. ■ 

Here, again, in order to apply the results of this section to formula (jlip we have 
to compute the values of the LT C{l,—B — 1,A), which can also be seen as 
a truncated gamma function. If B is an integer, C is easily computed by the 
binomial formula 

poo D| ^ \n 

C{l,-B-l,X) = l il + -fe-'''dx^^J2^-. 
If _S > is not an integer, ([77|) gives 

/oo P('R_|_1\ 1 

x^e-^-dx = ^^^p/ - iFi(l; 2 + B; A) 

r(^ + 1) 

XB + 1 

but then we have to rely on numerical analysis for computing the corresponding 
values of the confluent hypergeometric function 

iFAl: B + 2; X) = 1 + ^ 7 r 

' ' ^^{B + 2)...{B + n + l) 

and use (P71) . A good reference for such numerical consideration aspects can be 
found in Abad and Sesma (1995). 



r(B + 71 + 2) 

Tl— _ 



Practical conclusion for a — \,h — —\ — B . We fix a pair (c, d) such that 
either condition 1 or condition 2 of Proposition 12 holds. We choose a number A 
such that < Ai = A — c < A + c — A2. We compute numerically C(l; — 1 — i?, Ai) 
and C(l;— 1 — i?, A2). The mixing coefficient p is therefore determined by (fTTj) . 
The density on (0,oo) 

+ K(1;-1-B,A,) + ^1 - ^)C(1;-1-S,A2). 
has an increasing hazard rate. 
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